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Abstract 

Using a semiclassical continuum model of an electron in a deformable molecular 
crystal, some properties of multicomponent generalizations of the polaron- "vector 
polarons"- are elucidated. Analytical solutions for the case of two electronic bands 
coupled to two vibrational modes are given in detail. Within the model considered, 
the vector polaron can be classified by its wavefunction into several types and can 
have features that include: (1) a spatial variation in the electronic and vibrational 
character, and (2) low-energy internal degrees of freedom. For the case of electronic 
and vibrational degeneracy, local Jahn- Teller interactions can also lead to a novel 
spatiotemporal soliton, a long-lived excited state of the many-electron system sta- 
bilized by the conservation law resulting from degeneracy. 
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1 Introduction 



Many of the complex solids currently under intensive theoretical and exper- 
imental study have common attributes. High superconductors, fuUerene- 
based materials, and colossal magnetoresistive manganites all feature elec- 
tronic and vibrational degeneracies that are a consequence of the spatial sym- 
metries of the crystal structure. These systems typically have multiple elec- 
tronic and structural phases competing at low temperatures. And there is 
ample evidence [1] that the electron-phonon interaction and the resulting po- 
laron quasiparticles play an important role in shaping many of the physical 
properties of these materials. 

These observations serve as motivation for the theoretical study of the effects 
of generalized electron-phonon coupling in multiband systems. Using a semi- 
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classical continuum model of an electron in a deformable molecular crystal, 
the effects of electronic and vibrational degeneracy (and near degeneracy) on 
polaron structure have been investigated. With the inclusion of local Jahn- 
Teller couplings within the degenerate manifolds, a multicomponent polaron 
emerges as a solution of the semiclassical equations of motion [2] . 

Experimentally the richer vibrational and electronic structure of this "vector 
polaron" might well be probed using STM to image the charge density in the 
vicinity of a vector polaron trapped at a defect. This would be similar to the 
STM studies of the surface of an anisotropically gapped superconductor with 
impurities [3], where the order parameter symmetry has a great effect in shap- 
ing spatial features of the local electron density of states in the vicinity of an 
impurity, and this can be seen through measurements of the local conductivity. 



2 The Model 



Consider a molecular solid with two orbitals per molecule: l^/^i) and \4>2)- The 
molecule is also taken to have two vibrational modes: a breathing mode with 
amplitude 0o and a symmetry-breaking mode with amplitude 0i. The orbitals 
are taken to transform under elements of the molecular point group so that 
the symmetry-breaking mode couples to the difference in orbital densities, a 
local E ® j3 Jahn- Teller term. 

The continuum Hamiltonian for this system is taken to be 



H = Hel + Hph + HjT 

with 

H,i = I d'^x'i'^ (-^^oV + Acts) * 



(1) 

(2) 
(3) 



and 



HjT — 



(4) 



where \E'(a;) is a two-component field operator. The upper component destroys 
an electron from the first orbital at x, while the lower component destroys an 
electron from the second orbital at x. In the absence of distortions, the band 
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constructed from {ipi) states is offset in energy by 2 A from tlic band con- 
structed from 1-02) states. The limit A ^ corresponds to doubly degenerate 
bands. 

A continuum approximation to the local E ® (3 model studied by Hock et al. 
[4] is recovered in the limit (^o ~^ 0, while the continuum Holstein molecular 
crystal model [5] is obtained in the limit gi 0. Hence, the model of Eq. 1 
combines the Holstein model and the E ® j3 model. 



3 Bogoliubov-de Gennes Equations 



The Bogoliubov-de Gennes (BdG) equations for this model can be found by 
expanding the field operator in terms of a new real-space basis. Envelope 
functions t'n}^i that modulate the molecular orbitals are introduced, 
and the field operator becomes 



* = E !"^n (5) 



where c„ destroys an electron with amplitudes Un{x) and Vn{x) for the two 
molecular orbitals. A necessary condition (unitarity) on the amplitudes is then 

j d'^x{ulUn' + V^Vn') = ^nn' (6) 

The BdG equations for the envelope functions are then found to be 

(^-^V^ A - 5(000 - Un = enUn 

-^V'^ - A - go(l)o + gi(l)ij Vn ^ e„Vn (7) 



From energy minimization, the distortions must satisfy the following self- 
consistent equations 



00= — ^{\Unf + l^^nH 

</'i=-E(KP-kP) (8) 

where the sums are over all occupied states. 
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For a system with a single electron, only one term contributes to the sum. 
Eqs. 8 can be substituted into the BdG equations to obtain the two coupled 
nonlinear Schrodinger equations for the envelope functions u, v previously ob- 
tained in Ref. [2] 



(^-^V^ - u\u\'^ - r]\v\'^^ u= {e- A)u 

(^-^V^-?7|m|^-i/|v|^^ v = (e + A)v (9) 
where u = gl/ kq + gl/ and rj = gl/ Kq - gj/ Ki. 

The single electron might also be thought of as occupying a conduction band 
or partially occupying a set of conduction bands. Provided the polarization 
effects from the filled valence band are small, the "frozen" valence band ap- 
proximation can be made, and the identical BdG equations result. 

The spatial coordinates and the envelope functions can be scaled 



C = \/2\e-A\x 

(10) 



ri \ / u ( u 



T2/ Y |e - A| 
to yield the following equations 



(V^ - 1 + + I3rfj ri = 
(V' -uj^ + rl + (3rl) r2 = 



fill 



where o;^ = (|e| — A)/{\e\ + A) and (3 — r)/u. 



Eqs. 11 have been previously studied in many different contexts such as non- 
linear optics (propagation of light in birefringent optical fibers [6] ) and atomic 
physics (spinor Bose- Einstein condensation [7]) where the resulting Ginzburg- 
Landau equations have the same form. Although the mathematical study of 
systems of coupled nonlinear Schrodinger equations is still in its infancy, a 
set of known solutions in ID and 2D to these equations that describe the 
components of the vector polaron wavefunction will now be summarized. 
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4 Vector Polcirons in ID 



For one spatial dimension, three regimes can be identified where analytic solu- 
tions to Eq. 11 that are stable and localized are known to exist. In general, it is 
required that e < and |e| < A. 77 and v are restricted here to be non-negative. 
Thus, < /3 < 1. 



^.1 Manakov type 

The case of = 1 and u = 1 is special, as it corresponds to an integrable 
system. Manakov [8] studied this system of equations in another context. The 
symmetry of these equations gives rise to an infinite set of degenerate solutions 
that can be labeled by a continuous parameter < ^ < 27r. 

ri{^) = a/2 cos 6* sech ^ 



ui = 1 corresponds to exact degeneracy (A = 0), while l3 = 1 results from no 
coupling to the symmetry-breaking mode {gi — > 0). Consequently, the energy 
of the state is independent of orbital density differences. Thus the Manakov 
solution contains a zero-energy internal degree of freedom, corresponding to 
periodic changes in the local orbital mixing. 



4-2 Equal Amplitude type 

One situation that is exactly soluble is the "equal amplitude" case where 
ri(^) = r2{C,)- For stable, localized solutions, the parameters must be further 
restricted so that a; = 1. In this case, Eqs. 11 reduce to 



The solution to this scalar nonlinear Schrodinger equation is the well-known 
bell sohton 



^2(0 = V^sinO sech ^ 



(12) 




(13) 




(14) 
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The continuous symmetry present in the Manakov case is broken for /3 7^ 1. 
However, rotations of the orbital mixing still serve as a low-energy degree of 
freedom for {3 close to 1. 



4.3 Wave-Daughter Wave type 



For a; < 1, corresponding to non-degeneracy of the orbitals, it is possible 
to have one envelope function dominate throughout the core of the localized 
wavefunction. In this wave-daughter wave approximation, the amplitude of the 
daughter wave r2 is much smaller than that of the primary wave ri, leading 
to the following approximate set of equations 



3^ - 1 + 



-uj'' + (3rl]r2^Q (15) 




The top equation is recognized as the scalar nonlinear Schrodinger equation 
with a solution 

ri(^) = sech ^ (16) 



The primary wave provides an attractive potential for the daughter wave, and 
the Poschl-Teller equation for the daughter wave results 

+ 2/3sech2e) r2 = (17) 



This is known [9] to have a single bound state for < /3 < 1, with 

r2(C) = a sech'^ ^ (1^ 



and an eigenvalue constraint that uj = (\/l -|- 8(3 — l)/2. 

Within these cases amenable to analytical solutions, it is seen that (1) vector 

polaron states can support low-energy internal degrees of freedom, correspond- 
ing to pseudo-rotations of the electron amplitude in the orbital space; and (2) 
vector polaron states can have spatial variations in the orbital mixing and in 
the character of the distortion, as illustrated in the wave-daughter wave type. 
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5 Vector Polcirons in 2D 



Vortex vector solitary wave solutions are known to exist in 2D. Consider solu- 
tions to Eqs. 11 of the form ri(p, 0) = exp(imj0). The following coupled 
equations result 



where Ap = \Plp)' Vortex solutions to Eq. 19 may be labeled by the 

topological charges |mi,m2). 

The simplest case is the |0, 0) solution, where both envelope functions are ra- 
dially symmetric. Using an analytical variational method, the existence of a 
localized |0, 0) solution has been shown [10]. Unfortunately this solution has 
also been shown to be hnearly unstable, as are all of these vortex solitons. Re- 
lated multipole vector solitons [12] are also linearly unstable in this simplified 
model. 

However, similar solutions that are linearly stable have been found for slightly 
modified systems of equations that include either saturable nonlinearities [10] 
or nonlocal nonlinearities [11]. Saturable nonlinearities would physically corre- 
spond to the addition of proper vibrational anharmonicity of higher orders in 
this model, while nonlocal nonlinearity would correspond to gradient terms in 
the elastic energy. Consequently there are indications that semiclassical mul- 
ticomponent polaron states of the vortex and multipole varieties stably exist 
in generahzed models of 2D systems. 



6 Jahn-Teller Soliton 

Another related continuum model that employs a local Jahn-Teller coupling 
can also lead to a stable soliton, a state involving a lump of excess electronic 
charge trapped by a time-dependent pseudo- rotation of the elastic field [13]. 
These sohtons fall in the class of solutions known as non-topological Q-balls 
[14] that have been invoked in models of cosmological phase transitions in the 
early Universe [15]. 

Consider a continuum model with three electronic bands linearly coupled to 
two degenerate elastic modes. Take the case of where the electronic bands 




(19) 
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transforms as x and y and a trivial orbital singlet; the two elastic modes 
transform as x and y. The system may be described by the following Hamil- 
tonian 



H^He + Hph + HjT (20) 
where 



//e= E ^dM~^l'^'^m + W5^ml,l^l^|;J (21) 

m=0,±l '' ^ 

Hph^ E fd^x{\dt<P^\^ + U{\<Pm\') (22) 
5 ^ / d^x{(j)_^ijlil:o + H.c.) (23) 

™ LI ^ 



m=±l 



m=±l 

where W is the energy splitting between the doublet bands and the singlet, 
and m labels the states in the axial angular momentum basis. 

Such a system has a continuous symmetry whose origin lies in the elec- 
tronic and vibrational degeneracy. By virtue of Noether's theorem, a conserved 
charge Q results from this continuous symmetry provided it remains unbroken. 
The constraint of fixed Q provides stability to the soliton state. 

For large Q, a variational ansatz called the thin wall approximation [14] has 
been found to give good results in calculating the sohton energy. The ampli- 
tude of the elastic deformation is taken to be constant over a finite region 
of space with volume V, while vanishing outside that region. This deforma- 
tion itself, however, is not constant; rather, it undergoes pseudorotations in 
the two-fold degenerate space. It is straightforward to show that the Noether 
charge Q is linearly proportional to the pseudorotation frequency cu and to 
the square of the amplitude of distortion [13] 

Q ^ 1/ (24) 



Thus, for sufficiently large band splitting W, the energy of this thin- walled 
soliton is 



E^{Kph + U- Ejt)V 

^-uj-'ev +u{e)v -^^ev 



Q 



2 



2ev 



+ U{e)V (25) 
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where Kph and U are the kinetic and potential energy density of the time- 
dependent distortion. Ejt is the energy density from the electron-phonon 
interaction, no is the electron density, and U = U{C,^) — ^^p^C^- The last step 
follows from using Eq. 24 to eliminate cu in favor of Q and ^. 

As Coleman [14] showed, the resulting energy can be minimized subject to 
fixed Q to find the soliton energy and its size. 





2U 




Q 



(26) 

V = (27) 

Thus provided that the minimum value of 2U /^^ occurs at ^ ^ 0, this Jahn- 
Teller soliton cannot decay by emitting phonons that carry away the charge 
Q, since the energy per charge is lower for the sohton. This criterion sets the 
requirements on the amount and type of anharmonicity that are necessary to 
stabihze the soliton. 

In summary, it is noted that many of the solitons, both topological and non- 
topological, that are found in multicomponent Ginzburg-Landau theories can 
be used to describe vector polarons within the semiclassical theory outlined 
here. Experimental STM studies of these states trapped at defects might 
be used to reveal the non-trivial spatial structure of these polarons. Vector 
polarons, together with the many-electron Q-ball-like excited state, serve as 
examples that systems containing electronic and vibrational degeneracy can 
support novel states through local Jahn-Teller interactions. 
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